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Abstract 

A simple class of unitary renormalization group transformations that force hamiltonians towards 
a band-diagonal form produce few-body interactions in which low- and high-energy states are de- 
coupled, which can greatly simplify many-body calculations. One such transformation has been 
applied to phenomenological and effective field theory nucleon-nucleon interactions with success, 
but further progress requires consistent treatment of at least the three- nucleon interaction. In 
this paper we demonstrate in an extremely simple model how these renormalization group trans- 
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I. INTRODUCTION 



Wilsonian renormalization group transformations are designed to replace explicit coupling 
between disparate distance or energy scales with effective interactions in which disparate 
scales are decoupled [l|. In recent work we have employed a simple unitary renormalization 
group transformation to study the nucleon-nucleon (NN) interaction 0, y, HI]. This trans- 
formation is a simplified version of Wegner's flow equations 0] and one of a much larger 
class of Similarity Renormalization Group (SRG) transformations developed by Glazek and 
Wilson [f| . The transformation leads to NN potentials for which calculations of few-nucleon 
binding energies and other observables converge rapidly 0,1,0 . However, further progress 
will require the consistent treatment of at least the three-nucleon interaction. 

The basics of the unitary evolution are simply stated. Consider a unitary transformation 
U 8 of a hamiltonian H, 

H s = U s HUl (1) 

where s is a continuous flow parameter. We take s = for the initial value with Uq the 
identity transformation, so H is our input hamiltonian. We want to choose U s so that 
H s is diagonalized (band-diagonalized in more realistic cases) within a specified basis as 
s —>■ oo, which will realize the desired decoupling of low- and high-energy states. Direct 
differentiation shows that H s evolves according to 

= (2) 

with 

V(s) = = -r/t( s ) . (3) 

Choosing r](s) specifies the unitary transformation. A simple choice is 

r)(s) = [T,H s ], (4) 
where T is a fixed matrix (independent of s), which gives the flow equation, 

^ = [[T,H S ],H S }. (5) 

In Refs. 

as a 

T was chosen to be the kinetic energy; here we choose the hamiltonian for 
free particles in an infinite square well. For explicit calculations we employ a basis in which 
T is diagonal and it is in this representation that H s is driven towards band-diagonal form, 
as we will see below. 

A principal advantage of SRG transformations is that all operators are consistently trans- 
formed, which means that all observables are invariant. For the simple unitary transforma- 
tion this is obvious, simply because it is unitary. An additional advantage is that SRG 
transformations readily handle Fock space operators; indeed, Glazek and Wilson designed 
them to attack light-front quantum chromodynamics. Interactions that change particle 
number are not required for low-energy nuclear physics, but we do require the consistent 
evolution of all many-body operators. If we express the hamiltonian in terms of creation and 
destruction operators, it is evident that the commutators in Eq. ([5]) will generate many-body 
interactions even if H includes only a two-body interaction. In principle this could make 
practical calculations intractable, but in applications of interest (e.g., to nuclear physics) 
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we can choose transformations that maintain a hierarchy of many-body forces such that for 
sufficiently dilute many-body systems we only need to evolve few-body operators. 

In this paper we illustrate how a unitary SRG consistently evolves two-body and three- 
body interactions. We choose what may be the simplest possible example of such evolution, 
bosons in a two-level system in which only two states are coupled in each of the two-particle 
and three-particle sectors. As a specific realization, we use a one- dimensional infinite square 
well with the Fock space truncated to allow only the two lowest modes. Despite the simplicity 
of the model, it exhibits many of the basic issues to be confronted in evolving few-body 
interactions and illustrates a diagrammatic formulation of the SRG equations that carries 
over to more realistic problems (such as nuclei). 

In the severely truncated model space, the largest sub-matrix we will encounter is 2 x 2; 
so we start in Sect. [II] by reviewing the unitary transformation for the simple case of 2 x 2 
matrices. In Sect. HHI we explicitly show in our simple example that the unitary evolution of 
hamiltonians containing two-body interactions must produce three-body interactions. We 
derive the explicit 2x2 matrices that represent the flowing two-body interaction, V2 S , and 
the positive parity part of the flowing three-body interaction, V 3s . In Sect. [IV] we develop 
simple diagrammatic rules for the production of evolution equations governing any iV-body 
interaction, Vn, and show that these readily reproduce the equations for Vi and V3 in our 
simple model. The diagrammatic representation suggests several approximations that are 
tested. We summarize and make connections to nuclear physics and other problems of 
interest in Sect. IVl 



II. 2x2 MATRIX EVOLUTION 

Let us start by considering a two-state system, with the hamiltonian represented by a 
2x2 hermitian matrix. In realistic many-body quantum mechanics problems, we frequently 
encounter a hamiltonian that naturally divides, H = T + V, with T a one-body operator. 1 
In this case, the unitary transformation 

H s = U S TU\ =T + V S (6) 

defines the evolved interaction V s . In terms of basis states where T \i) = E% \i), the flow 
equation is 

k 

Here the renormalization flow parameter, s, has the dimension of inverse energy squared. 
The first term on the right-hand-side drives V, and therefore H, towards band-diagonal form 
as s increases from zero. 

In our two-state example we will also take H = T + V and use eigenstates of T as a basis, 
with T\l) = Ei\l) and T|2) = _£?2 1 2 ) . In this representation, 



There is considerable freedom in the choice of T, it need not even be a one-body operator. V can be 
initially chosen to be a two-body operator, but the commutators will then automatically generate three- 
body operators, four-body operators, etc. 
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and the full hamiltonian is 



_ E x + V n {s) V 12 (s) \ 
s ~\ V 21 (s) E 2 + V 22 (s)J ' W 

where Vij(s) = (i|V^|j), which we take to be real. Any real symmetric 2x2 matrix can be 
expressed using the identity matrix, J, and the Pauli spin matrices, a x and a z , 

T = 81 + Sla z , (10) 

V s = cl + u z (s)a z + u x (s)a x , (11) 

where £ = \(E X + E 2 ), Q = \(E X - E 2 ), c = \(V 1X + V 22 ), u z (s) = \(V 1X - V 22 ) and 
Wx{s) — V12 = V 2X . Since / commutes with all other matrices, £ and c drop out of the 
flow equations and c is therefore independent of s. So, we need only consider the Pauli 
spin matrices a z and a x to understand the real 2x2 hamiltonian. If an additional angle is 
introduced using a y , it turns out to be another flow constant that produces no qualitative 
change in the way the hamiltonian flows. 

Dropping the constants £ and c, which can simply be added back in after the flow 
equations are solved, we work with 

V s = w z {s)a z + u x {s)a x , (12) 

and 

T = Qa z . (13) 
With the above matrices, the unitary renormalization group flow equation (jSj) becomes 

-^[u z (s)a z +LO x (s)a x ] = 40[w^(s)<7 a - (SI + u 2 (s))u x (s)a x ] , (14) 
or, projecting out coefficients of the Pauli matrices, 

^z(s) = 4&k%{s) , J-^ 8 ) = ~ 4fi ( fi + uj z (s))to x (s) . (15) 

These equations illustrate a tremendous advantage of such unitary transformations in that 
the exact renormalization group equations are only second order in V. Note that if >> 
|a; z (s)|, oo x (s) is obviously driven to zero exponentially. 

The equations are easily solved an alytically We in troduce the angle 9 and set SI + u z = 
uj cos 6 and u x = u sin 8, where to = ^ (SI + oo z ) 2 + uj\. With this change of variables we find 
that uj is a flow constant, so only the angle 9 depends on s. The flow equation becomes: 

d6 

— = -Attojsm6. (16) 
as 

Integrating this equation we obtain: 

tan = tan (M) e — . (17) 

Thus 9 is driven exponentially either to zero or to 7T, depending on the sign of SI. In either 
case the matrix is driven to diagonal form, with the states in the same order as in Sla z . 
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III. BOSONS IN A SQUARE WELL TRUNCATED TO TWO MODES 



In this section we consider a simple concrete example that lets us examine SRG evolution 
in two-boson and three-boson sectors of Fock space. We choose T to be the hamiltonian for 
non-interacting bosons in an infinite square well where —L/2 < x < L/2, 



We work in the T eigenstate basis, 



T=^ + V well . (18) 
2m 



2 . fmv(x-L/2)\ 

sin 1 - 1 } , (19) 



L \ L ) 



and truncate the boson Fock space to allow only the ground and first excited state modes. 
We will see that this drastic truncation reduces the two- and three-boson problems to 2 x 2 
matrix problems, as solved in Sect. [Ill 

After truncating to two modes, the complete set of eigenvalues of T are 

We can build many-body interactions using Fock-space operators and below we develop a 
Fock space diagrammatic analysis exploiting the simple algebra of these operators, but for 
pedagogical reasons we first compute matrix elements of the hamiltonian with two-body 
and three-body interactions mechanically. The interactions are given by their matrix ele- 
ments between unsymmetrized two- and three-boson states, which must be related to matrix 
elements between properly symmetrized states. We choose an interaction that preserves par- 
ity and need only consider the three-dimensional two-boson space and the two-dimensional 
positive parity three-boson space to display the method. 

It is convenient to use several bases for the many-boson problem, which we distinguish 
by the type of bracket used. Unsymmetrized three-boson states, for example, are \n1n2n3), 
in which rij's denote the individual particle square-well states. We also use symmetrized 
states \NiN 2 ), in which Ni is the number of bosons in the n = 1 state and N 2 the number 
in the n = 2 state. The only two-boson symmetrized states are 1 20) = \a), 1 02) = \b), and 
1 11} = |c) and we only need the positive parity symmetric three-boson states, 1 30) = \a) 
and 1 12) = \(3). 

We display the two-boson and three-boson states we need in Tables UHlVt listing their 
free energy and parity. We also list the symmetric states' representations in terms of unsym- 
metrized states, which we need to compute the initial hamiltonian. If the initial two-body 
interaction conserves parity, parity will be conserved in the evolution, so even and odd parity 
states will not mix. In general, any operator that commutes with T and H will commute 
with H s and result in an explicit symmetry. Parity conservation simplifies both the two- 
and three-boson problems, reducing them to 2 x 2 matrix problems. 

We will see that we must allow for an explicit three-body interaction, V3, without which 
the transformation cannot be unitary. If we set V3 = at s = 0, we will see that the 
unitary transformation generates a three-body interaction that maintains unitarity in the 
three-boson sector of Fock space. That is, we need V s = V% a + Vz s here. In realistic cases we 
need to choose the transformation so that the induced part of these many-body interactions 
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TABLE I: Two-boson unsymmetrized states 



\nin 2 ) 




parity 


|H) 


2Ei 


+ 


|22) 


8E 1 


+ 


|12) 


5E 1 




|21) 


5Ek 





TABLE II: Two-boson symmetrized states 



\NiN 2 ) 


label 


k 1 i n 2)-basis 


En u n 2 


parity 


|20) 


\a) 


|H) 


2Ei 


+ 


|02) 


\b) 


|22) 


8E 1 


+ 


|H) 


\c) 


73[|12) + |21)] 


5E 1 





TABLE III: Positive parity three-boson unsymmetrized states 



\ni n 2 n 3 ) 


E 


parity 


|1H) 


3E 1 


+ 


1 122) 


9E 1 


+ 


|212) 


9E 1 


+ 


|221) 


9E 1 


+ 



TABLE IV: Three-boson symmetrized states 



\N ± N 2 ) 


label 


b 1 i n 2"'3)-basis 




parity 


|30) 


\a) 


|111) 


3E 1 


+ 


|12) 


\P) 


J-[|122) + |212) + |221)] 


9E 1 


+ 


|21) 


It) 


^[|112) + |121) + |211)] 


6Ei 




|03) 


\S) 


|222) 


Y1E X 





TABLE V: Initial matrix elements of V 2 

(a\V 2 \a) = (b\V 2 \b) =~3g~ 
{a\V 2 \b) = 2g 
(c\V 2 \c) = 4g 
{a\V 2 \a) = 9g 
((3\V 2 \f3) = llg 
(a\V 2 \f3) = 2VSg 
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remains short-ranged, which for fermions should allow us to truncate the tower of many-body 
operators and accurately solve many-fermion problems with evolved few-body operators only. 
The calculation here is easily extended to the four-boson sector of Fock space, where 
first appears, and so on. Restricting our attention to the two-boson and three-boson sectors, 
we must solve 

^ = [[T, V 2s + V 3s ],T + V 2s + V 3s ] . (21) 

Solving this in the two-boson sector completely determines V 2s and the three-boson sector 
then determines V 3s . Since particle number is conserved, we simply need to take two-boson 
or three-boson matrix elements of the flow equation and insert a complete set of states on the 
right-hand-side to obtain coupled nonlinear equations for each independent matrix element. 
Two-boson matrix elements of the flow equation determine V 2s , 

^- s (a\V 2s \a) = -12E 1 (a\V 2a \b)(b\V 2a \a) , (22) 

^ s (b\V 2s \b) = 12£ 1 <&|y 2a |a)<a|V 2 ,|&) , (23) 

-r- s {a\V 2s \b) = -36E*{a\V 2s \b) + 6E X [(a\V 2s \a) - (b\V 2s \b)] (a\V 2s \b) . (24) 

There is no coupling between the positive and negative parity states for the interaction we 
will study and (c|V2 S |c) does not evolve because there is only one odd-parity two-boson state. 
We have seen in Sect. [Ill that these equations have an analytic solution; to apply it we need 
only remove the part of H that is proportional to an identity matrix and determine Q, lv z (0) 
and 0^(0) from what remains. 

The three-boson evolution equations are readily computed using eq. (12~TT) . 

^(a\V 2s + V 3s \a) = -\2E x {a\V 2s + V 3s \f3){f5\V 2s + V 3s \a) , (25) 

^(P\V 2s + V 3s \P) = 12£!</3|y 2s + V 3s \a)(a\V 2s + V 3s \(3) , (26) 

^-(a\V 2s + V 3s \/3) = -WEl{a\V 2s + V 3s \P) 



ds 



+ 6E 1 (a\V 2s + V 3s \a)(a\V 2s + V 3s \(3) 

-6E 1 (a\V 2s + V 3s \P)(P\V 2s + V 3s \P) . (27) 



Here again we have a 2 x 2 matrix problem, so it also can be solved as shown in Sec. [TT1 
One way to find V3 is to solve for V 2 + V 3 and then subtract matrix elements of V 2 that 
result from solving for V 2 in the two-boson sector. Instead, we will discuss a seemingly more 
complicated procedure that is required in more realistic examples where matrix elements of 
dV 2 /ds in the three-boson sector produce momentum delta-functions due to non- interacting 
spectators. 

We need symmetrized three-boson matrix elements of V 2s in terms of its symmetrized two- 
boson matrix elements, in order to embed the solution of eqs. (j2~2~]) -([2"I ]) in eqs. (1251) - (127j) . 
These are computed using unsymmetrized states given in Table IIII1 

V 2s \(3) = (b\V 2s \b) + 2(c\V 2s \c) , (28) 
V 2s \a) = V3(b\V 2s \a) , (29) 
a|^ 2s |a) = 3(a\V 2s \a) . (30) 
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We can now rewrite all of the three-boson matrix elements of V 2s in terms of its two-boson 
matrix elements, then use the solution to eqs. (I22l) - (l24l) to rewrite the three-boson equations 
for V 3s : 

^(a\V 3s \a) = -UE^ia^sl^^lV^a} + V3(a\V 3s \(3}(b\V 2s \a) 

+ (a\V 3 ,\P)(j3\V 3 ,\a)] , (31) 

^(P\V3 S \/3) = 12£ 1 [2(&|V 2 »(a|F 2s |&) 

+ V3(b\V 2s \a}(a\V 3s \(3) + V3(/3\V 3s \a) (a\V 2s \b) 

+ (P\V 3s \a)(a\V 3s \P)] , (32) 



^-W\V 3s \P) = -36El(a\V 38 \P) + 12E 1 V3((a\V 2s \a) - (c\V 2s \c)){a\V 2s \b) 
as 

+6E 1 [((a\V 3s \a) - (f3\V 3s \(3))V3(a\V 2s \b) + 

(S(a\V 2s \a) - (b\V 2s \b) -2(c\V 2s \c))(a\V 3s \(3)] 
- QE^V^a) - ({3\V 3s \(3)){a\V 3s \(3) . (33) 

There is no obvious simplification here, quite the contrary. Replacing matrix elements of V 2 
with their solutions from the two-boson sector complicates the equation, but we will return 
to this issue after developing a diagrammatic analysis. It is no longer obvious that this is 
another 2x2 matrix example and it is far more tedious to find an analytic solution, but 
once again the entire hamiltonian is driven to diagonal form, which requires V 2 and V 3 to be 
separately driven to diagonal form. 

For an explicit demonstration problem that can be exactly solved and readily extended, 
we choose a zero-range initial interaction: 

V 2 (xi,Xj) = 2Lg5(xi - Xj) . (34) 

Note that for g > this is repulsive. Matrix elements of V 2 between unsymmetrized states 
are: 

(nin 2 \V 2 \ 

HZ^a) — 9 [^0,rai— n2+ri3— 714 ,rti-n 2 -n3+n4 

"0,ni— n2+ri3+»i4 ^0,n\— n^— n^— 714 < ^0,ni+ji2+"3 — "4 
^0,rti+ri2— nz+rii $0,ni+n2— n^— 714] • (35) 

From these we readily compute the initial interaction matrix elements for symmetrized 
states, which are given in Table [Vj 

In figs. [T] and [2] we plot two-boson matrix elements of V 2s and three-boson matrix elements 
of V 3s , for weak and strong couplings. We contrast flow for g = —0.1 and g = —3.0, showing 
unitary renormalization group evolution as a function of s. The energy scale is set by 
choosing E\ = \. 

The trace of V 2 in the two-boson sector is fixed, so we always see one diagonal matrix 
element increase and the other decrease. Since the off-diagonal matrix element must die 
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FIG. 1: Hamiltonian matrix elements of V2 as a function of the flow parameter s in the two-boson 
sector for weak and strong couplings. 




10 4 10 3 10 2 10 1 10 4 10 3 10 2 10 



FIG. 2: Hamiltonian matrix elements of V3 as a function of the flow parameter s in the three-boson 
sector for weak and strong couplings. 



exponentially, it is not surprising to find similarity in fig. [T] of how V2 flows for strong and 
weak coupling. For g = —0.1 the rate of exponential convergence is controlled by the linear 
term on the right-hand-side of eq. [2H while for g = —3.0 the second-order term contributes 
significantly. This will be detailed in the next section. To compare this flow with that 
of a 2 x 2 matrix, note that here Q = —3Ei and for weak coupling the exponent will be 
approximately 36s, which is why we see convergence by the time s ~ 0.1. Convergence 
improves as the coupling increases, but the range of s over which convergence is seen will 
not change drastically until g is orders of magnitude larger than E\. 

V3 displays more interesting behavior. Its trace in the three-boson sector is not fixed, 
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only the trace of V2 + V 3 is fixed. The trace of V2 in the three-boson sector varies because 
of its evolution in the two-boson sector, so the trace of V 3 varies to restore unitarity. Here 
we see (/3\V 3s \j3) becoming much greater than (a|T^ a |a) for both weak and strong coupling, 
because {{3\V 3s \P) is fed by a term second-order in V 2 , which we will identify with a tree 
diagram in the next section. 

For g = -0.1, ((3\V 3s \(3) grows to 0(g 2 ) while (a\V 3s \a) only grows to 0(g 3 ). (a\V 3s \0) 
is initially comparable to (/3|T^ s |/3) because it is also fed by an 0(V 2 2 ) tree-level interaction, 
but it is then driven exponentially to zero with all off-diagonal matrix elements. Additional 
features of this flow will be discussed in the next section. 

-10 r ' J 



3 

> . 
a 

so 

T3 ' 
O 
X3 
I 

<u 

(H - 

3 




-35 - 



FIG. 3: The two lowest eigenvalues of the three-boson systems as a function of s with the full 
interaction and when the three-body interaction is omitted. 



One measure of the importance of V 3s is the relative error in the three-boson eigenvalues 
when V 3s is not included; the errors arise because V2 S alone can not produce unitary flow 
in the three-boson sector. The eigenvalues for these two cases using g = —3.0 are shown in 
fig. [3J The eigenvalues are constant when the full interaction is included ("It's unitary!"), 
while omitting V 3s leads to a small error for the ground state and a 50% error for the excited 
state. These errors are consistent with the missing matrix elements, (a|V 3s |a) and (l3\V 3s \f3), 
shown in fig. [2J For g = —0.1 the errors are too small to be easily visible. 

Of course, there is no reason to restrict the initial hamiltonian to include two-body 
interactions alone. The unitary flow of V2 S is not altered if we add an initial three-body 
interaction, but that of can be altered drastically. None of the flow equations derived 
above change, only the initial values of V3 matrix elements change. Off-diagonal matrix 
elements are still driven to zero and diagonal matrix elements yield the correct eigenvalues. 

We illustrate how a three-body force added from the start affects unitary renormalization 
group evolution using another zero-range initial interaction: 

V 3 (xi, Xj, x k ) = (2L) 2 g 3 5(xi - x s ) 5(xj - x k ) . (36) 

The initial matrix elements we need are (a|V3 S |a) = I0g 3 , (/3\V 3s \fl) = 6g 3 and (a|V3 S |/3) = 
5#3. Note the relatively large prefactors. 
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10 4 10 3 10 2 10 



s 

FIG. 4: Hamiltonian matrix elements of V3 as a function of the flow parameter s in the three-boson 
sector for weak coupling, g = —0.1 and g% = 0.01. 

In fig. H]we display typical perturbative evolution, choosing g = —0.1 and g 3 = 0.01 so 
that Vi and V3 have comparable perturbative effects. The dominant effect of the perturbative 
three-body force is to both shift and split the three-boson energies. The off-diagonal matrix 
element decays exponentially to zero with little additional effect on the energies. 




FIG. 5: Hamiltonian matrix elements of V3 as a function of the flow parameter s in the three-boson 
sector for two sets of strong couplings. 



In fig. Owe display two examples of non-perturbative evolution, both with g = —3.0 as 
considered above, with g 3 — 0.1 an example where V2 dominates the evolution and g% = 1.0 
an example where V3 dominates in the three-boson sector. These should be contrasted with 
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FIG. 6: Hamiltonian matrix elements of V3 as a function of the flow parameter s in the three-boson 
sector for two additional sets of strong couplings. 

fig. [2] in which g 3 = 0. Perhaps the most interesting feature to emerge as increases is 
that V3 splits the three-boson energies with a different sign for intermediate values of #3. 
Figure [6] shows that as #3 is further increased, the order of its three-boson splittings reverts 
to the perturbative ordering. This simple model actually displays a remarkably broad range 
of interesting, analytically understandable behavior as a function of the two couplings and 
our examples are not exhaustive. 

IV. FOCK SPACE DIAGRAMMATIC ANALYSIS 

While it is straightforward to generate coupled evolution equations for many-body in- 
teractions by computing explicit matrix elements of the operator evolution equations, this 
requires inserting a complete set of iV-boson states to determine and does not immedi- 
ately expose disconnected interactions in which spectators do not participate. If one seeks a 
numerically tractable differential equation for the evolution of Vn in momentum representa- 
tion, these disconnected matrix elements are best avoided because they contain momentum 
delta-functions for the spectator states. One of the advantages of a diagrammatic analysis is 
that it isolates these matrix elements as disconnected diagrams, which are exactly cancelled 
as long as V^r-i is fixed by solving the SRG equations in sectors with fewer than N particles. 

We want to develop simple diagrammatic representations of [[T, V], T + V], to generate 
explicit unitary flow equations for matrix elements of V. We must choose which many- 
boson matrix elements to represent. The matrix elements that immediately lead to simple 
diagrammatic rules are symmetrized, but not normalized for states containing more than 
one boson in a given mode. Since we use the same states on both sides of the unitary flow 
equation, normalization does not matter until we compute the explicit hamiltonian. 

Eigenstates of T define particle creation and annihilation operators, |n) = a^|0), with 
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[a m , ajj] = 5 mn , in terms of which, 

\i,j,k, ...} = a\a]a\--- |0) , (37) 

where i, j, k, . . . label external and internal legs in the diagrams. We associate an eigenvalue 
of T, with any line labelled by i. For two modes, 

T = E\a\ai + E 2 a\a 2 , (38) 



and for any number of modes, 



T = ^E i a\a i . (39) 



V2 is then specified by {zj'lV^lfc/} and V 3 is specified by {ijk\V 3 \lmn}: 

1 \ 2 



r2= (^) ^24al{ij\V 2 \kl} ai a k , (40) 
V 3 = (^j) ^2 4 a Wk^ k \ V 3\ lmn } a na m ai ■ (41) 

ijklmn 

In our simple example above, these indices took on only the values 1 and 2, but this constraint 
does not simplify the diagrammatic rules and we need impose it only at the end to verify 
that the rules produce the same equations as derived in the previous section. 

We are now in a position to develop a complete set of diagrammatic rules for the right- 
hand-side of the flow equations for {ij|V2|fc/} and {ijk\ V 3 \lmn}. In fig. Owe show a diagram- 
matic representation for the basic elements of these flow equations. All terms are built from 
V N ,V N = [T,V N ] tmdV N = [[T,V N ),T], where {ij\V 2 \kl} = (Ei + Ej — E^ — Ei){ij\V 2 \kl}, 

and similar relations are easily found for all matrix elements of and Vjy. External legs 
on these diagrams are labeled with symmetrized labels, so we will not count diagrams as 
"topologically distinct" if they result from a permutation of indices associated with the legs 
of a single vertex. The redundancies from these permutations cancel the (1/2!) 2 in V 2 and 
the (1/3!) 2 in V3, except in diagrams with loops, where an iV-loop diagram is weighted by 
l/(iV + 1)!. Diagrams with one loop are assigned a factor of 1/2, etc. 




FIG. 7: Labelled vertices represent matrix elements of Vn, Vn = [T, Vn], and Vjy = [[T,Vn],T]. 

In fig. [H]we show a schematic representation of the diagrams that appear in the SRG flow 
equations in the two-boson and three-boson sectors. Diagrammatic rules are easily deduced 
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ds 
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r/.s 



FIG. 8: Complete diagrammatic representation of dV2/ds and a schematic representation of the 
diagrams for dV 3 /ds. For the latter, there are eighteen tree diagrams, six one-loop diagrams, and 
two two- loop diagrams. 

from the equations and are readily generalized to any sector. For clarity we expand eq. (1211) : 

^f + ^T = [[T,V 2s ],T] + [[T,V 3s ],T] 

+ [[T, V 2s ),V 2s ) + [[T, V 2s ), V 3s ] + [[T, V 3s ],V 2s ] + [[T, V 3s ],V 3s ] 
= % + F 3 + { ([F 2 , V 2 ] + [V 2 , V 3 ] + [V 3 , V 2 ] + [V 3 , 7 3 ]) 

-(V 2 ^V 27 V 3 ^V 3 ]) \ . (42) 



Diagrams represent symmetrized matrix elements of these equations. Every line carries a 
label and products of V and V result in contractions that equate indices on lines between 
the two vertices. Once indices are assigned to all lines, vertices represent matrix elements of 
V 2 , V 3 , V2 and V 3 , with indices matched to lines. There is a sum over all internal indices, 
so even tree diagrams represent sums. One- loop diagrams appear in the evolution of V 2 , 
while both one- and two-loop diagrams govern the evolution of V 3 . In general, tree diagrams 
through (N — l)-loop diagrams appear in the evolution of Vn- 

The structure of eq. (14"2"|) guarantees that V 2s is completely determined by two-boson 
matrix elements, with V 3s then determined by three-boson matrix elements. Comparing 
eq. ( 1421) with fig. El note that V 2 and V 3 appear together in the equation but there are 
separate diagrammatic equations for their evolution. Two-boson matrix elements of eq. fj4"2j) 
produce the diagrammatic equation for dV 2 /ds, while three-boson matrix elements produce 
the diagrammatic equation for dV 3 /ds. When three-boson matrix elements of eq. (j4"2j) are 
computed, the matrix elements of dV 2 /ds on the left-hand-side cancel against diagrams on 
the right-hand-side that are identical to the two-boson diagrams with a spectator line that is 
not connected to either vertex. As long as dV 2 /ds solves the two-boson evolution equations, 
these disconnected diagrams can be dropped on the right-hand-side of the diagrammatic 
equation and only dV 3 /ds appears on the left-hand-side. 

There are two reasons that disconnected diagrams never appear. Diagrams in which 
both vertices are not connected to one another by at least one line vanish because of the 
commutator structure, [V, V}. If there are no contractions between creation and annihilation 
operators in these two vertices, the two terms from the commutator produce identical matrix 
elements that exactly cancel. 
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As discussed above, diagrams in which an external line is not connected to any vertex, 
such as three-boson diagrams in which a single spectator line appears with the one-loop 
interactions between the other two bosons as displayed below in fig. [9j simply produce 
copies of the flow equations from lower sectors of Fock space. So, if V 2 satisfies the two-boson 
flow equations, disconnected three-boson matrix elements of dV^/ds cancel the disconnected 
matrix elements of [Va, V 2 \. We isolate the flow equation for dV^/ds simply by dropping 
these disconnected diagrams and setting dV 3 /ds equal to the sum of connected diagrams, 

including the diagram for V3. 

So, V2 is determined by two-boson matrix elements, because V3, V4, etc. will not appear 
in this sector and all diagrams that do appear in the two-boson sector also appear in the 
three-boson sector with one disconnected line, in the four-boson sector with two disconnected 
lines, etc. Again, if V2 satisfies the flow equation in the two-boson sector, no disconnected 
diagrams survive in the flow equation for V3 in the three-boson sector. In the three-boson 
sector, there are tree diagrams involving two powers of V 2 , one-loop diagrams involving one 
power of V2 and one power of V3, and two-loop diagrams involving two powers of V 3 . V 2 is 
already determined by the two-boson evolution equation and V3 is the only new interaction 
appearing in the three-boson evolution equations, so the three-boson evolution equations 
determine V3. 

This type of analysis is easily generalized to higher sectors because of the simple second- 
order structure of the exact unitary flow equation. V4 is determined by the four-boson 
evolution equations. These include tree-diagrams with one power or V2 and one of V3, one- 
loop diagrams with one power of V2 and one power of V4 or with two powers of V3, two- loop 
diagrams with one power of V3 and one power of V4, and three- loop diagrams that are 
second-order in V4. 




FIG. 9: Diagrammatic representation of eq. (03]). 



In fig. [9] the complete set of labelled two-boson diagrams is shown, from which we obtain 
-j-{ij\V 2 \kl} = {ij\V2\kl} + ^J2({^\^2\mn}{mn\V 2 \kl}~{tj\V2\mn}{mn\V2\kl}) . (43) 

mn 

Substituting for V2 and V 2 we find 

^-{ij\V 2 \kl} = -(E id - E kl ) 2 {ij\V 2 \kl} + ^ J2( E v + Ekl ~ 2E mn ){ij\V 2 \mn}{mn\V 2 \kl} , 

mn 

(44) 

where Eij = E,i + Ej. It is straightforward to show that when sums are restricted to the 
lowest two modes, eqs. (T2^1) - (f24l) are reproduced. 

In fig. [TU] the complete set of labelled three-boson diagrams is shown. This equation is 
for the evolution of V3, with V 2 solving the two-boson equations above and no disconnected 
diagrams appearing as a result. There are eighteen distinct tree-level diagrams, twelve one- 
loop diagrams and two two-loop diagrams. We do not explicitly display all possible external 
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+ 8 i, j <-> k and m,n «-> I permutations 




+ 2 m, n <-> / permutations 




FIG. 10: Diagrammatic representation of eq. (|45p . The notation i, j «-> A; means i «-> j or j <-> A;. 

legs in the diagrams or the equations but list them as permutations. The resulting equation 
is 

-7"{*j'A;| V 3 |Zmn} = {zj/cl V^l/mri} 
as 

+ / ] ({i j '[ V2 1 fo} {pfe | V2 1 wn} — |14|^}{p^|^2|wn} 
+ 8 z, j «-» fc, m, n <-> I permutations) 
+^^2({ij\V 2 \pq}{pqk\V 3 \lmn} - {ij|F 2 |pg}{pgA;|V3|Zmn} 

pq 

+ 2 i, j «-» permutations) 
+^5Z({^'fc|V r 3 |/pg}{pg|V r 2 |mn} - {ii/c|y 3 |/pg}{pg|y 2 |wn} 

+ 2 m, n <-> Z permutations) 
+— ^^({zjfe|y3|pgr}{pgr|V3|Zmn} — {zjfc|V3|pgr}{pgr|K 3 |lmn}) . (45) 

pgr 

Most of the terms on the right-hand-sides of these equations are given by permutations of 
the indices in the terms explicitly displayed. For example, in the tree diagrams in fig. ITOl 
one particle line on the left is connected to a different vertex than the other two. We 
must explicitly add each permutation of the index assigned to this line, obtaining additional 
diagrams in which this index is i instead of k and j instead of k. Permutations of the two 
lines connected to the same vertex are considered "topologically identical" to the diagram 
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shown, so they produce no new diagrams. This reflects the fact that {ij '\V2\pq} = {ji\V2\pq}. 

In front of each sum over N indices there is a factor 1/N\, so to each iV-loop diagram 
we associate this factor of I/AH. These sums are unrestricted; if we replace them with sums 
over distinct permutations the factor of 1/JV! drops out. 

It is again a straightforward but tedious exercise to show that truncation to the two 
lowest modes yields eqs. (l3"Ti) - (13"3l . To match equations, note that \a) = ■s/l/2 |11}, \b) = 
v/T/2 1 22} and |c) = 1 12} are the two-boson states we need, while \a) = a/1/6 1 1 1 1 } and 
= a/1/2 1 122} are the only positive parity three-boson states we need. 

If the hamiltonian has no three-body force at the start of SRG evolution, one will be pro- 
duced. Initially, the tree diagrams, which are second-order in V 2 , dominate. Next V3 will 
suppress far off-diagonal matrix elements while terms of order V2V3 will be more important 
near the diagonal, with the two- loop diagram emerging last since it is order V^. Naturally, 
if V3 is not zero at the start of the renormalization group evolution, the relative magnitude 
of these contributions can change. The diagrams provide a tool for analyzing various contri- 
butions to the flow equations and seeking approximations, which might be essential in more 
realistic calculations. 



250 1 




FIG. 11: Matrix elements that contribute to ^(a|V2 S |6) in the two-boson sector for weak and 
strong coupling. (a|p2 S |6) is dominated by V2, with the second-order one-loop correction barely 
visible when g = —0.1. 

In figs. [TTHT5] we contrast various contributions to the evolution of V% and V3 for weak 
coupling, g = —0.1, and strong coupling, g = —3.0. Figure [TT1 displays contributions to the 
off-diagonal matrix element (a|V2 S |&). In both cases the initial evolution is dominated by the 

double-commutator, V2, which always drives off-diagonal matrix elements to zero exponen- 
tially. For weak coupling the one-loop, 0(V 2 2 ), contribution remains negligible throughout 

the evolution, while for strong coupling it becomes comparable to that of V 2 and helps drive 
(a|V2 S |6) to zero. This qualitative behavior is what is seen when this transformation is ap- 
plied to realistic nucleon-nucleon interactions [1, but it is possible to find cases where 
the one-loop contribution first opposes V 2 and forces (a|V2 S |&) to grow before V 2 finally takes 
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FIG. 12: Matrix elements that contribute to ^{o\V2 S \o) in the two-boson sector for weak and 
strong coupling. Only the one-loop diagram contributes. 
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FIG. 13: Matrix elements that contribute to -f^{(y\V^ s \(3) in the three-boson sector for weak and 
strong coupling. Tree diagrams initially dominate the evolution of (a|V3 S |/3), allowing V3 to grow 
and then drive (a|V3 S |/5) to zero. Here again, the one-loop and two-loop contributions are barely 
visible for g = —0.1 and the two-loop contribution remains negligible even for strong coupling. 
This can change if V3 is not zero initially. 



over [9|. This happens when the running cutoff encounters a bound state threshold and 
we make no attempt to discuss such interesting exceptions in this paper, although they are 
inevitable in three-body systems that display the Efimov effect (lfj, 111 . 

Figure [121 shows the sole one- loop contribution to the evolution of the diagonal matrix ele- 
ment (a I V2s|a). Since the trace of H s is conserved in each sector of Fock space, contributions 
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FIG. 14: Matrix elements that contribute to ^(a|V3 S |a) in the three-boson sector for weak and 
strong coupling. No tree diagrams contribute to (a|V3 S |a), and the relatively small one-loop con- 
tributions dominate over the two-loop contributions in both cases. 
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FIG. 15: Matrix elements that contribute to ^(/?|V3 S |/3) in the three-boson sector for weak and 
strong coupling. Tree diagrams dominate the initial evolution of (/?|V3 S |/?), with relatively small 
one-loop contributions and negligible two-loop contributions in both cases. 



to (b\V2s\b) simply have the opposite sign. The only change from weak to strong coupling is 
the magnitude of the one-loop contribution. This simplicity is partially a reflection of the 
drastic truncation of Fock space but a similar exponential fall-off in strength is guaranteed 
in more realistic examples as long as the transformation drives the matrix to band-diagonal 
form. Not only do the far off-diagonal matrix elements flow to zero exponentially, the diago- 
nal matrix elements evolve rapidly at first and then slow as they approach their asymptotic 
values. Again, there are exceptions to this simple evolution if bound-state thresholds are 
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FIG. 16: Relative errors in the three-boson eigenvalues for weak and strong coupling when V3 is 
simply dropped. 



encountered (Sj. 

In figs. [T3HT51 we turn to terms that contribute to the evolution of V3. In fig. [13] we 
immediately see that tree diagrams, which are (9(V 2 2 ), dominate the initial evolution of 
(a|V3 S |/3) for both weak and strong coupling. This is guaranteed when V3 itself starts at 
zero, but the subsequent evolution differs for weak and strong coupling. For weak coupling, 

once (a|V3 S |/3) itself has any appreciable strength, V3 builds and drives this off-diagonal 
matrix element back to zero. One-loop and two-loop contributions remain negligible for 
all s when the coupling is weak, simply because these terms are suppressed by powers of 
g. At strong coupling it is the one-loop contribution that builds to shut off the growth of 
(a|V3 S |/3) and then drive it to zero, while V3 and the two-loop contribution remain small but 
non-negligible near the end of the evolution. 

The evolution of (a|V3 S |a) and (/3|V3 S |/3) are not simply connected by a trace constraint 
because V2 also contributes to the trace of H s in the three-boson sector. There is no tree 
diagram contributing to (a|V3 S |a), so it remains small in comparison to {(3\V^ s \j3) for both 
weak and strong coupling. We see in fig. [TH that the one-loop diagram dominates for 
both weak and strong coupling, but this contribution changes sign during the evolution 
of (a|V3 S |a) for strong coupling. Figure [TB] shows that the tree diagram dominates this 
evolution, with other contributions remaining negligible throughout the evolution at weak 
coupling and only the one-loop contribution becoming significant at strong coupling. 

In fig. [16] we show the fractional error in the three-boson eigenvalues E a and Ep when V3 
is simply dropped. This error remains smaller than 0(g 2 ) for weak coupling, which is the 
order of the tree diagrams that are ignored, and grows to 50% for strong coupling. In both 
cases it is Ep that is most sensitive to V3. 

In fig. [T7] we show these same fractional errors after approximating V3 by evolving with 

only the tree diagrams and V3. This significantly reduces the errors in Ep from simply 
dropping V3 for both weak and strong coupling. This is consistent with the observations 
above about the small relative strength of the one-loop and two-loop contributions. Finally, 
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FIG. 17: Relative errors in the three-boson eigenvalues for weak and coupling when V3 and tree 
diagrams are used to evolve V3. 




FIG. 18: Relative errors in the three-boson eigenvalues for weak and strong coupling when tree 
and one- loop diagrams are used to evolve V3. 



in fig. [18] we show these errors when V3 is evolved using both tree and one-loop diagrams, 
dropping only the two-loop diagrams. This further reduces the errors for weak coupling 
but not for strong coupling, where we have seen in fig. [13] that the two-loop contribution, 
although smaller than the one-loop contribution, becomes comparable to the contribution 

of V3 throughout the evolution. 

It should be clear that this analysis is easily extended to the four-boson sector, where 
V4 must be added. For weak coupling, we should again find that tree diagrams which are 
^(^2^3) will dominate in the evolution of V4. Since there are no diagrams contributing to 
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V4 built only from two powers of V2, V4 is fed entirely by induced interactions and might 
remain small even for strong two-body interactions, but explicit calculations are required to 
find what happens at strong coupling. 



V. SUMMARY 

In this paper, we have used a simple two-level model system to illustrate how a unitary 
SRG consistently evolves two-body and three-body interactions. We have developed a di- 
grammatic treatment of the SRG equations that applies generally and have shown how it 
can be used with a variety of examples in our model. However, such a model might seem 
too schematic to offer insight into realistic problems. In fact, the steps to the physics prob- 
lem that originally motivated this work, that of describing atomic nuclei from microscopic 
interactions, are not so large. 

The low-energy nuclear many-body problem has been one of the most important problems 
in physics for over seventy years, but it has stubbornly resisted a complete solution. In recent 
years, however, it has become clear that a Wilsonian renormalization group perspective leads 
to new conceptual and practical ways to attack this problem. The primary tool Wilson 
introduced to renormalization group formalism is a renormalization transformation that 
lowers a well-defined hamiltonian resolution, a tool that is ideally suited to finding universal 
interactions. Feed any interaction constrained by low-energy data into a well-designed RG 
transformation, it produces a universal low-energy interaction that is decoupled from the 
ambiguous (and irrelevant) high-energy components of the input interaction. 

The same RG transformation that Wilson employed in his first non-perturbative RG 
calculation [12j has recently been applied to a wide variety of modern nucleon-nucleon inter- 
actions 13|, ll4| and every one of them was transformed into a nearly universal low-energy 
effective nucleon-nucleon interaction, called V\ ow k- It turns out that Vj ow fc also drastically 
reduces the complexity of the nuclear many-body problem but requires the consistent evo- 
lution of three- nucleon interactions [l6| . 

Realistic nucleon-nucleon (NN) interactions include strong short-range repulsion and a 
strong tensor force. These lead to far off-diagonal strength in the momentum representation 
of the nucleon-nucleon interaction, so that large matrices are required to represent this 
interaction. The size of these matrices increases rapidly with particle number. Vi ow fc requires 
drastically fewer matrix elements to accurately reproduce all low-energy NN scattering data, 
and this reduction exponentiates as we move to the many-nucleon problem [3]. In contrast 
to conventional NN interactions, which are highly nonperturbative, V\ ow k leads to nearly 
converged nuclear matter calculations at second order in many-body perturbation theory. 
The three-nucleon interaction is an essential ingredient in saturating nuclear matter close to 



the empirical density and binding energy [15|. Even more, the problem is not well-posed until 
consistent few-body interactions are found, as the nucleon-nucleon interaction on its own 
does not maintain unitarity in even the three-nucleon problem and violations are comparable 
to nuclear level spacings [7|. 

The SRG transformation we employ produces basically the same universal nucleon- 
nucleon interaction as the Vj ow fc evolution and shares the same positive features. Perhaps the 
most important advantage of the type of SRG transformation used here is that it "automat- 
ically" allows us to consistently evolve two- and few-nucleon interactions. Just as we have 
seen in our model, if the 3-body force is not included in the evolution, few-body observables 
are cutoff dependent and errors can be large. This article represents a first step towards the 
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computation of a fully consistent low-momentum three-nucleon interaction. In addition to 
the three-nucleon interaction that is required by unitarity when the two-nucleon interaction 
flows, there is an intrinsic three-nucleon interaction that must be added; current versions 
include various phenomenological and chiral effective field theory potentials. Ideally the 
SRG transformation will again reveal a universal low-momentum, three-nucleon interaction. 

In principle, four-nucleon interactions can also be included and may be needed for quan- 
titative nuclear structure, particularly for heavy nuclei. The three-loop diagrams in the 
SRG equations represent nine-dimensional integrals and the spin-isospin algebra also grows 
in complexity. But we have seen in our model that even for strong couplings we can some- 
times accurately approximate the SRG evolution equations by dropping such multi-loop 
contributions and we may even obtain better than 10% accuracy using tree diagrams alone; 
10% accuracy for an already small correction might be sufficient. There are no integrals in 
these tree diagrams, so we need "only" deal with the spin-isospin algebra for four nucleons. 
But this is idle speculation until the three-body interaction is computed accurately. With 
a three-nucleon interaction in hand we can directly test the important of V4 by studying 
violations of unitarity in 4 He, for example. 

The basic procedure for evolving two-nucleon and three-nucleon interactions with a uni- 
tary RG transformation is the same as that illustrated in this article. The square well must 
be removed, bosons must be replaced by nucleons, one dimension goes to three dimensions 
and the initial hamiltonian must include realistic inter-nucleon interactions. Removing the 
square well and replacing bosons with fermions is uncomplicated, while moving up to three 
dimensions and using realistic interactions requires more work but no major conceptual 
developments. 

It is important to realize that in nuclear calculations, we do not intend to use the trans- 
formation to diagonalize the hamiltonian as we have in this article. In realistic many-body 
problems we run the flow parameter s only to the point where low- and high-energy degrees 
of freedom are decoupled in the few-nucleon sectors. Running the transformation beyond 
this point would produce long-range few-body forces and strong many-body forces; we would 
effectively be forced to solve both few- and many-nucleon problems completely. We want 
to simplify the few-body forces that are then employed in many-nucleon calculations, not 
diagonalize the hamiltonian in every sector of Fock space and produce the large tower of 
many-body forces that would be required to do this in the many-body sectors. When we 
achieve sufficient decoupling in the two-nucleon sector, we compute the three-nucleon in- 
teraction using methods illustrated in this article and then switch to conventional few- or 
many-body methods using this partially diagonalized hamiltonian. See Refs. ^ and for 
details. 

Perhaps the most interesting immediate extension of the model we present is to remove 
the square well and let T simply be the kinetic energy. Instead of discrete sums, we obtain 
integrals over continuous momenta, but these can be approximated using a discrete mesh, 
so we once again have sums over discrete indices. iV bosons interacting via a two-body 
delta-function potential in one dimension produce one bound state whose wave function 
and energy are known analytically [13]. Using momentum eigenstates we must compute the 
evolution of {ki, ki \ V2 \ k[, k' 2 } and {ki, k^, k% | V3 | k[, k' 2 , k' 3 }. These can then be used in 
a numerical basis function calculation that should converge increasingly rapidly with basis 
size to the analytic result as the unitary transformation is run. With this discrete basis, the 
evolution equations are almost identical to those in this article. 

The basic unitary renormalization group equation we employ in this article, eq. (jSj), is 
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the same in three dimensions as in one. Replacing bosons with fermions means only minor 
changes in the diagrammatic rules for the evolution equations, with the familiar minus signs 
from fermion exchange and no symmetry factors associated with loops. Exploiting angular 
momentum algebra and choosing appropriate Jacobi coordinates for three-nucleon states 
introduce technical complications, but these appear in a form almost identical to what is 



faced when solving three-nucleon Faddeev equations [18| . The SRG equations governing the 



evolution of the nuclear three-body force are only a scaled-up version (larger matrices) of our 
model problem, as opposed to the multi-channel complications of the Faddeev equations. 

In this article we have chosen one simple operator for T, but there are many possible 
choices and this freedom can be exploited to improve convergence in many-body calculations. 
The calculations in this article are easily modified for any T that can be defined by diagonal 
matrix elements in the square well basis. T can itself depend on the flow parameter. Wegner 
jH| advocated using the diagonal matrix elements of H s itself, which corresponds in the 2x2 
matrix model to choosing T = u z (s)a z . Wegner's transformation is guaranteed to drive 
the momentum representation of the hamiltonian to diagonal form, while other choices can 
produce divergent matrix elements and stall. These considerations are discussed in the 
context of simple models in Ref. [9|]. 
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